Revisiting the Dolinar Receiver through 
Multiple— Copy State Discrimination Theory 
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We consider the problem of discriminating between two quantum coherent states by interpreting 
a single state like being a collection of several successive copies of weaker coherent states. By means 
of recent results on multiple-copy state discrimination, it is possible to give a reinterpretation of the 
Dolinar receiver, and carry out a quite straightforward analysis of its behavior. We also propose and 
investigate a suboptimal detection scheme derived from the Dolinar's architecture, which is shown 
to slightly outperform some other near-optimal schemes available in literature. 
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I. INTRODUCTION 

Discrimination between two non-orthogonal quantum 
states is a fundamental issue in quantum mechanics and, 
in particular, in quantum communications. From a the- 
oretical point of view, the problem was completely ana- 
lyzed and solved by Helstrom [lj, which found the opti- 
mal measurement operators, and the corresponding cor- 
rect detection probability (Helstrom bound), for both 
pure and mixed quantum states. Unfortunately, also 
for pure states, often the optimal measurements do not 
correspond to quantum observables that are easily mea- 
surable, so that the experimental implementation of the 
optimal discrimination is a very difficult task. 

For the case of two coherent states of a traveling single 
radiation mode, in 1973 Dolinar [2j proposed an adaptive 
measurement scheme, based on a combination of photon 
counting and feedback control, that precisely achieves the 
Helstrom bound (see also Q)- However, since the scheme 
requires a very precise control of an optical-electrical 
loop, only recently the Dolinar's idea has obtained a sat- 
isfactory practical implementation [J|. 

Recent years have seen an increasing interest for adap- 
tive measurements from both a theoretical and an experi- 
mental point of view, also for optical phase measurements 
and estimation (see Q and references therein). A no- 
tably interesting theoretical result has been obtained by 
Acin et al. [6| for discrimination between pure quantum 
states, when multiple identical copies of a quantum state 
are available. They proved (see also p}) that in this case 
the Helstrom bound can be achieved by local adaptive 
measurements applied to single copies. The result is par- 
ticularly attractive in that it offers a useful insight into 
the Dolinar's approach for discrimination between coher- 
ent states. In this paper we discuss the strict connection, 
already recognized in |5|, between the ideas underlying 
measurements of multiple copies of pure states [6| and 
the Dolinar receiver 0|. 

Acin et al. [6j considered the discrimination between 
pure quantum states I70) and I71), when n identical 



copies of an unknown state are given. Formally, the prob- 
lem consists in discriminating between the pure states 
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in the tensorial product Hilbert space < H® n , where IK 
is the Hilbert space spanned by the single copies I70) 
and I71). Of course, also in this case one could apply 
the Helstrom theory and find an optimal collective mea- 
surement in J-C® n achieving the Helstrom bound. Unfor- 
tunately, collective measurements are difficult to realize 
experimentally. With the adaptive approach suggested 
in |6j, the experimenter performs on each copy a local 
measurement which is optimized on the basis of the re- 
sults of the measurements on the previous copies. The 
surprising enough conclusion is that the optimized local 
measurements achieve the Helstrom bound, exactly the 
same as the optimal collective measurement. 

The discrimination between two coherent states of a 
single-mode harmonic oscillator (without loss of gener- 
ality 1 7) and I —7)) presents a difficulty similar to that 
of collective measurements on multiple copies. Namely, 
the Helstrom theory gives optimum measurement vectors 
that are linear superposition of I7) and | — 7) and do not 
correspond to any measurable observable. On the other 
hand, owing to their peculiar properties, the coherent 
states 1 7) and | — 7) of duration T can be thought as se- 
quences of shorter and weaker modes of duration T/n, 
namely, 



l7> = 
-7> = 



(2) 



Moreover, as n increases, and the average number of 
photons per copy goes to zero, the optimal Helstrom 
measurement on each copy may be conveniently approx- 
imated by a displacement followed by a photon detec- 
tion. Then, in principle, we may think to apply the 



multiple-copy adaptive measurement to the segmented 
quantum states $2$. As n goes to infinity, it appears 
natural the transition to the Dolinar scheme [2J], with a 
continuous time-varying displacement controlled by the 
photon counting results. 

This paper is organized as follows. In Section HI] 
we illustrate the optimal multiple-copy measurement 
paradigm from a novel point of view leading in a natural 
way to the continuous-time extension as in the Dolinar 
receiver. In Section Mil the feasibility of near-optimal 
discrimination of weak coherent states is discussed. In 
Section lTvl we extend the optimal multiple-copy adaptive 
measurements to coherent states. In particular, we derive 
the theory of the Dolinar receiver in a simple way, avoid- 
ing the cumbersome machinery of dynamic programming. 
Finally, in Section [V] we propose a suboptimal simplified 
version of the Dolinar's scheme. 



II. MULTIPLE COPY ADAPTIVE 
MEASUREMENT 

Measurement strategies for discrimination between 
multiple copies of two pure quantum states were dis- 
cussed by Acin et al. [6fl (see also p}). Alice, accordingly 
to the binary input symbol a G {0, 1}, chooses between 
two pure states I70) and I71) in the Hilbert space "K with 
probability qq and q± = 1 — qo, respectively. (Without 
loss of generality we assume q$ > qi). Then, Alice sends 
Bob n identical copies of the chosen state, corresponding 
to the states ([T]), which are pure states in the tensorial 
product Hilbert space !K®™. Bob performs a measure- 
ment on the system, trying to guess the original state 
with maximum correct detection probability. 

The optimal result is given by the well-known Hel- 
strom bound [l| 
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where 



A=|(ao|ai)| = |(7o|7i)r=X ? 



(4) 

and x = I (To Iti) I is the overlap coefficient of the sin- 
gle copies. Bob may achieve this bound by a global 
measurement using suitable von Neumann projectors 
Bo = |/?o)(A)| an d Bi = |/3i)(/3i| over the product space 
!K® n . Unfortunately, the optimum measurement vectors 
|/3o) and \j3\) turn out to be non separable linear superpo- 
sition of the pure states \<xq) and |ai), i.e., an entangling 
measurement which is hard to implement experimentally. 
The problem of optimizing local adaptive measure- 
ments had been tackled by Acin et al. [6| and it may 
be formalized in the following way. Let us assume, with- 
out loss of generality, that the single copies are given by 



I70) = cos#|a;) + sin#| 
J7i) = cos#|a;) — sin0| 



(5) 



where \x) and \y) form an orthonormal basis of the 
Hilbert space "K spanned by the states I70) and I71) and 
the overlap coefficient is given by 



X = (7o|7i) = cos 29 



(6) 



Assume that the local measurement orthonormal vectors 
on the A:-th copy, k = 1, . . . , n, 



|/jfco) = cos <f>k\x) + sin<^ fc |y) 
|/jfci) =sin<£ fe |a;) -cos<j> k \y) 



(7) 



are completely specified by the measurement angle <f>k- 
Let Zk € {0, 1} be the outcome of the fc-th measurement, 
which is also assumed as the result of a provisional de- 
cision. The adaptive optimization problem consists in 
finding a starting measurement angle <\>\ and a recursive 
rule 



<f>k = fk(zi, ■ ■ ■ ,Zk-i) 



(8) 



in such a way that the final outcome z n gives the correct 
detection with maximum probability. In the recursion 
([8|) the information gained by the previous k — 1 mea- 
surements is exploited in order to optimize the choice of 
the next measurement angle. 

This appears to be a dynamic programming problem 
[8| and, as such, it had been dealt with and solved in [6fl . 
The main results, surprisingly simple, are summarized as 
follows: i) the optimal solution of the dynamic program- 
ming approach gives correct detection probability coin- 
ciding with the Helstrom bound ((3]), so that the global 
optimal measurement may be replaced by more easily 
implementable local measurements; ii) the problem re- 
duces to a bayesian updating problem (see also p]]) with 
recursive relation fa = fk(Zk-i) so that the new optimal 
measurement angle depends only on the outcome of the 
last measurement; iii) if the (k — l)-th result is Zk-i = i, 
the optimal measurement angle <pk is the solution of the 
Helstrom optimization problem obtained replacing the a 
priori probabilities q and q\ with the a posteriori prob- 
abilities P [a = 0| Zk-\ = i] and P [a = 1 1 Zk-i = i], re- 
spectively. 

In particular, after the measurement on the (fc — l)-th 
copy, the provisional correct detection probability coin- 
cides with the Helstrom bound on k — 1 copies, namely, 
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The next measurement angle <f>k for the k-th copy is 
chosen to maximize the probability of correct detection 
under the assumption that the a priori probabilities q 
and qi are replaced by the corresponding a posteriori 
probabilities of the input simbol, given the last result 

Zk-i- These turn out to be P [a = i\zk-i = i] — Pc , 
i = 0,1. Finally, the measurement angles are given by 



bk = — arctan 



1 
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(10) 



if Zfc-i = and 7r/2 — <pk if Zk-i — 1. A simple proof of 
these results is given in the Appendix. 1 

The optimum local adaptive measurement can be sum- 
marized by the following step-by-step procedure. 

1. From the overlap coefficient \ an d the input prob- 
abilities go and qi compute the two sequences of 
measurement angles 
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2. Start with angle <f>\ if qg > 1/2 (and tt/2 — (pi oth- 
erwise). 

3. Use the angles of the first sequence (|11[) as long as 
the measurement result is 0. 

4. Change angle sequence every time the result 
changes and accept z n as the global result. 

We will show in the sequel that this paradigm is mim- 
icked in a continuous time version by the Dolinar receiver. 

As a further comment, we note that the multiple-copy 
optimization requires the discrimination between two hy- 
potheses in the 2™-dimensional Hilbert space 3-f® n , so 
that the optimal solution is not uniquely defined. On the 
contrary, the Helstrom solution discriminates between 
the hypotheses in the restricted subspace 'Kq spanned 
by |ao) and \u\). Of course, each optimal measurement 
in 3i® n , once projected in J{ , returns the Helstrom pro- 
jectors. In particular, in the above procedure, to each se- 
quence z\ , . . . , z n of results it corresponds a measurement 
vector \fii) <£) . . .<g) \fi n ) in "K® n with measurement angles 
chosen in the sequences (fTTj) . It can be easily verified that 
the 2™ measurement vectors in "K® n are orthonormal and 
they globally give a von Neumann projective measure. 



III. COHERENT SINGLE-COPY 
MEASUREMENT 

Now we consider the possibility of applying the 
multiple-copy adaptive approach to the discrimination 
between two coherent states segmented like in ((5). Let 
us suppose that Alice prepares a single copy of binary 
coherent states. Without loss of generality, we can as- 
sume \ao) = I7) and \a\) = | — 7), with 7 real, as in the 
Binary Phase Shift Keying (BPSK) modulation scheme, 
with overlap coefficient 



X = |(a |ai)| 






(12) 



1 Note that the above results cannot be extended to multiple copies 
of mixed states Q . 



where j 2 represents the average number of photons in 
each state. A straightforward application of the Helstrom 
theory leads to the Helstrom bound 
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but the corresponding measurement vectors are difficult 
to implement. Then, several near-optimal, but simpler 
to implement, detection schemes had been devised in the 
past. 

A well-known solution is given by the Kennedy receiver 
[10J . Bob applies a displacement D{— 7) to the Alice state 
and tests the resulting state with a photon counter. The 
displaced states become \a' ) = |0) and |a^) = | — 27). 
The photon counting detection corresponds to measure- 
ment projectors Ho = |0)(0| and Hi = I — Hq. Decision 
a = 1 is accepted if the photon counter clicks, otherwise 
a = is chosen. The correct detection probability reads 



Pk =qo(a' \U a \a' ) +<7i(ai|niK) 



(14) 



An improved version [ll| of the Kennedy receiver em- 
ploys a displacement D(—/3) to be optimized, so that the 
correct detection probability becomes 



Pik =qo(l - /3|n | 7 - 0) + qi (7 + /3|Hi | 7 
=c? e- ( ^ )2 + gi (l - e-^ 2 



(15) 



By nulling the derivative with respect to /3, we find that 
the displacement quantity /3q maximizing (|15p satisfies 
the transcendental equation 
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(16) 



that can be numerically solved, and the corresponding 
value of Pik evaluated. 

In Fig. [T] the performance of the Kennedy receiver and 
of the improved Kennedy receiver are compared with the 
Helstrom bound. (The figure also includes the simplified 
Dolinar receiver that will be introduced in Section [Vj . 
For large values of 7 the improvement obtained by op- 
timizing the displacement /3 is negligible. On the other 
hand, as 7 goes to 0, both the Helstrom bound and Pik 
approach similar values. Indeed, as shown in the fig- 
ure, the performance of the improved Kennedy receiver 
strictly approximates the Helstrom bound for weak co- 
herent states. 



IV. A SIMPLE APPROACH TO THE DOLINAR 
RECEIVER 

The above considerations suggest that, for n large 
enough, such that copies of weak enough coherent states 
are obtained, the optimum measurements on the seg- 
mented states ^ can be well approximated by suitable 
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Figure 1: Symbol error probability in the case of equiproba- 
ble states (qo — qi = 1/2). The performance of the consid- 
ered detection strategies are reported and compared with the 
Helstrom bound, which is achievable by the optimal Dolinar 
receiver. Emphasis is given to the quantum limited region 
corresponding to very weak coherent states. 



displacements and photon counting. In other words, the 
sequences of measurement angles (|11|) may be reinter- 
preted as sequences of displacements. Then, the transi- 
tion to the continuous time-scheme depicted in Fig. [5J 
appears to be natural. 

The input field i^(i), < t < T, corresponding to the 
coherent state | ± 7) , is represented by 



ip(t) = ±ipe 



i27rfot 



(17) 



where /o is the optical frequency and T is the pulse du- 
ration. The mean number of photons arriving to the 
detector is given by 
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\iP(t)\ 2 dt = i) 2 T 



(18) 



The detector subtracts from the input field a time- 
varying field generated by a local laser with envelope 
chosen between either u (t) or t*i(£), accordingly to the 
value of z(t), a binary signal with possible values and 
1, giving the provisional decision at time t. By mimick- 
ing the behavior of the optimal multiple-copy detection, 
we assume that the decision signal z(t) changes at any 
photon arrival at the counter. Then, the optical signal 
at the photon-counter has envelope either ±ip — uq (t) or 
±tp — ui(t), depending on the value of z{t). Moreover 
z(T) is assumed to be the final decision. 

The mathematical problem is to choose the functions 
Uo(t) and u±(t) that maximize the correct detection prob- 
ability 



P c = P[z(T) 



(19) 



The problem can be solved by means of standard photon 
counting statistics. 
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Phot, on 
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Figure 2: Block diagram of the Dolinar receiver: The received 
signal with envelope ±i/> is displaced by a quantity —U z (t), 
where z(t) £ {0, 1} is the temporary estimation at time t. 
The value of z(t) alternately changes from to 1 at each single 
photon detection of the displaced signal ±i[>—u z ( t y The shape 
of the feedback signal u z i t \ is also changed accordingly. 



Let us assume that a — 0, so that ip(t) = -0e i27r ^°*. 
Then, t he p rocess z(t) can be interpreted as a telegraph 
process [12J alternately driven by non-homogeneous Pois- 
son processes with rates 



X(t) = |^-uo(i)| 5 



and 



H(t) = |^- Wl (t)| 2 . (20) 



A simple application of the properties of Poisson pro- 
cesses gives the evolution of the conditional correct de- 
tection probability po(t) = P [z(t) = 0\a = 0]. In fact, 
let N(t, t + At) denote the number of photon arrivals at 
the counter in the interval (t,t + At], therefore 

Po (t + At) = 

= P[z(t) =0,N(t,t + At) =0|o = 0] 

+ P [ z(t) = 1, N(t, t + At) = 1 1 a = ] + o(Ai) 
= P[N(t.,t + At) =0\z(t) = 0} Po (t) 

+ P [N(t, t + At) = 1 1 z(t) = 1 ] [1 - po(*)] + o(At) 
= [ 1 - \(t)At] Po (t) + n{t)At[l - po(*)] + o(Ai) . 

(21) 



Hence, the differential equation 
S Po (t) 



Po(*) 



St 



li(t) - [X(t) + n(t)]po(t) 



(22) 



follows. In a similar way for pi(t) = P [z(t) = l\a = 1 
we get 



p , 1 (t) = 6 -^- = m-[Kt)+Ht)}pi(t) 



with 
\(t) 
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(23) 



and 



m 
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u (t)\ 2 . 

(24) 

If our search is confined to symmetric solutions, 

namely, Ui(t) = — Uo(t) we get X(t) = X(t) and fl(t) = 

/J,(t), and the correct detection probability satisfies the 

differential equation 



Kit) 



SPcjt) 

St 
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= QoPo(t) +qip[(t) 
[A(i)+/x(i)]P c (t). 



(25) 



On the basis of the results on multiple- copy measure- 
ments we expect that, for some choice u{t) of the enve- 
lope of the feedback signal Uo(t), the provisional correct 



detection probability P c {t) is exactly equal to the Hel- 
strom bound applied to the interval (0,t), namely, 



1 



Pc{t) = - [l + Vl - 4g 9ie- 
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(26) 



By substituting the above expression in (|25[) , and defined 
R(t) = s/l - Aq« qi e-^ 2t , we get 



v 



R(t) 

= *P 2 + u 2 (t) + 2^u{t) - [V> 2 + u\t)] [1 + R(t)} 

(27) 



and after some algebra 
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(28) 



coinciding indeed with the Dolinar's solution (see 
also QJ]). 



V. A SUBOPTIMAL RECEIVER 

The Dolinar receiver requires a time-varying feedback 
signal uo(t), whereas both the Kennedy receiver [lj| and 
its improved version [ll| make use of a constant fixed dis- 
placement (3 leading to a much simpler implementation. 
Therefore, it is worthwhile to consider a simplified ver- 
sion of the Dolinar receiver where the feedback signal is 
constrained to have a constant fixed envelope Uo(t) = /3. 
Such a setting would mean that only phase modulation, 
and specifically phase inversion, is required, whereas the 
optimal receiver has to also employ an amplitude modu- 
lator capable of generating an optical signal with shape 
defined by (|28|) . which decays with t but it is divergent 
about t = for the important case of equiprobable states 
qo = qi = 1/2 . 

Hence, by substitution in (|24|) . i.e., setting A(t) = A = 
\ip — /3| 2 and fi(t) = \x = |V> + /?| 2 , with the initial condi- 
tion P c (0) = (ft)) we get the following final probability of 
correct decision 



c[ j_ 2% 2 +/3 2 



qo ~2 
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(29) 

The optimized value of j3 can be found by numerically 
solving the following transcendental equation, which is 
obtained by nulling the derivative of P C (T) made with 
respect to j3 

/?T(V, 2 + P 2 ) [(2a, - 1)(V' 2 + (3 2 ) - 2W] e-^+^ T = 
= ^ 2 -/? 2 )sinh((?A 2 +/3 2 )T). 

(30) 

From Figure [T] we note that the simplified Dolinar re- 
ceiver slightly outperforms the improved Kennedy re- 
ceiver. Therefore, its experimental demonstration, and 
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Figure 3: Intensity |/3| 2 of the fixed displacement D(—/3) for 
the considered receivers; setting go = 9i = 1/2, T — 1. 



a study of its robustness to the presence of possible im- 
pairments, is an interesting task for future contributions. 
In Figure [3] the intensity of the displacements for dif- 
ferent schemes is reported. It can be noted that, as con- 
firmed by Figure 1, it is for very weak coherent signals 
that the simplified Dolinar receiver and the improved 
Kennedy receiver are particularly attractive, since, with 
the increasing of the signal strength, they both performs 
a displacement similar to the one applied by the orig- 
inal Kennedy's proposal, and both (fTSj) and (|29|) ap- 
proach (H3J). 



VI. CONCLUSIONS 

A whole single coherent quantum state can be inter- 
preted as a succession of many (possibly infinite) weaker 
copies of the same state, so that a different interpretation 
of the quantum discrimination task can be given. Such 
a view may provide further insights to well-established 
problems and give rise to novel detection solutions. 

We have proposed an analysis of the behavior of the 
Dolinar receiver based on recent findings in the field of 
multiple-copy state discrimination. With such an ap- 
proach it has been possible to provide a quite intuitive 
explanation, and simple mathematical derivation, of the 
Dolinar receiver. We also proposed a suboptimal simpli- 
fied detection scheme that employs a photon counter, a 
phase inverter and, contrary to the Dolinar's solution, a 
constant-envelope displacement. 
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Appendix: Optimality of the multiple copy 
adaptive measurement 



The proof is done by induction. Using measurement 
vectors ([7]) with angle <f>\ , the correct detection probabil- 
ity after the first measurement turns out to be 



^ (1) -9o|(^io|7o)| 2 + gi|(Mii|7i>| 2 

= g cos 2 (# - 4>i) +gisin 2 (# + (f>i) 



(A.l) 



and, as it can be easily verified |6j] , it is maximized by <\>\ 
given by (jTDJ) for k = 1, with provisional correct detection 
probability given by 



f 



pP = 2 [i + \A - Mmr 



(A.2) 



Then, the result is proven for k = 1. In particular, 
simple computations give the a posteriori probabilities 



P[a = i\zi=i] =P { c l) . 

Now, suppose that the result holds true for fe— 1. From 
the inductive hypothesis, the provisional correct detec- 
tion coincides with the Helstrom bound, the adaptive 
measurement up to the (k — 1)— th copy coincides with the 
optimal global measurement and the a posteriori proba- 
bilities are 



P[o 



• Uk-l 



P?" 



(A.3) 



If these probabilities replace go an d <?i in the expression 
(fTOf of the angle <pk and in (JA.2I) , one gets 



p (fe) = 



1 -AP. 



(fe-i) 



1-P 



(fc-i) 



(A.4) 



1 + V 7 ! - MoqiX 2k 
and the proof is complete. 
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